February 5, 2013 



m 
o 

(N 



Two Particle-Hole Excitations in Charged Current Quasielastic Antineutrino-Nucleus 

Scattering 

J. Nieves,"'^ I. Ruiz Simo,^ and M. J. Vicente Vacas'^ 

^ Instituto de Fisica Corpuscular (IFIC), Centra Mixta Universidad de Valencia- C SIC, 
Institutas de Investigacion de Patema, E-4-6071 Valencia, Spain 
^ Dipartimenta di Fisica, Universita di Trenta, 1-38123 Trenta, Italy 
Departamenta de Fisica Teorica and IFIC, Centra Mixta Universidad de Valencia- C SIC, 
Institutas de Investigacion de Patema, E-46071 Valencia, Spain 

We evaluate the quasielastic and multinucleon contributions to the antineutrino nucleus scattering 
cross section and compare our results with the recent MiniBooNE data. We use a local Fermi gas 
model that includes RPA correlations and gets the multinucleon part from a systematic many body 
expansion of the W boson selfenergy in the nuclear medium. The same model had been quite 
successful for the neutrino cross section and contains no new parameters. We have also analysed 
the relevance of 2p2h events for the antineutrino energy reconstruction. 
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Neutrino (antineutrino) charged-current (CC) scattering on nuclei without emitted pions is a crucial detection 
channel for many neutrino experiments ^1-3] and has been extensively studied by several theoretical groups [4- 2 b]. 
i-S^ , It has been customary to take for granted that most of those eventS"'^ could be attributed to the quasielastic (QE) 
scattering of the weak probe on a nucleon. Hence, the initial neutrino(antineutrino) energy could be approximately 
O |. determined from the energy and angle of the final lepton assuming QE kinematics. 

D ' This assumption soon led to some puzzling results. In particular, neutrino MiniBooNE data Q were quite chal- 
^ ' lenging showing substantially larger cross sections than the theoretical predictions [13, H^l • A simple way out was to 
modify in the calculations some standard hadronic parameters, such as the nucleon axial mass Ma to fit 

data. However, the large obtained values (e.g. Ma — 1.35 ± 0.17 GeV in Q) were in clear confiict with other electron 
I> . and even neutrino experiments which result in a value very close to one {Ma = 1.014 ± 0.014 GeV in (2^). 
CO ■ A more natural solution came from the incorporation of some standard nuclear effects such as RPA and multinucleon 
mechanisms. The importance of these effects in the non mesonic neutrino nucleus scattering was first explored in 
' Refs. [1, Q and later in Refs. [H, [l^, [l^. Some of these more complete models were found to describe well even the 
. I MiniBooNE double differential cross section while using standard values for all parameters of the calculations [l^ . In 
. these models, a substantial part of the cross section measured in Ref. [1] corresponds to events in which at least two 
^5 ' nucleons are emitted. 

The relevance of the multinucleon mechanisms has some unwanted consequences. Obviously, the neutrino energy 
reconstruction, based on the QE kinematics is not so reliable [1^ [25l - l27j and that implies larger systematic uncer- 
tainties in the experiments analysis. Furthermore, a clear separation of quasielastic and multinucleon events is not 
straightforward from the experimental point of view, even if final nucleons were detected, due to the energy thresholds 
; ^ , and to the strong final state interaction of the nucleons in their way out of the nucleus. It is evident the need for 
■ comprehensive and theoretically well founded models that provide a good description for single and multinucleon 
processes. Only with their help, one could approach the tasks of disentangling the genuine QE events and properly 
reconstruct the neutrino energy. Of particular interest are the atomic number and the neutrino energy dependence of 
the reactions, the quest for kinematical regions with minimal multinucleon contamination and the relative importance 
of multinucleon mechanisms for neutrino versus antineutrino induced reactions. In fact, the latter observable could 
be an stringent test for the available models [16i. i28j. Furthermore, if the contribution of multinucleon mechanisms is 
substantially different in neutrinos and antineutrinos, as predicted for instance in Refs. [^[l^ and this is not properly 
understood, it could lead to an asymmetry between v and P which could be misinterpreted as a consequence of CP 
violation. 

Our approach for the study of QE and multinucleon processes induced by neutrinos was developed in Refs. |i,[l[ll- 
The model starts from a local Fermi gas (LEG) picture for the nucleus. The purely QE contribution for both 
neutrinos and antineutrinos was studied in Refs. incorporating among other nuclear effects Pauli blocking, Fermi 
motion and the long range correlations (RPA). A correct energy balance is imposed using the experimental Q values. 
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^ After subtraction of some background, such as pion production followed by absorption. 
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Of the various options discussed in those papers, we wiU use here the fully relativistic model of Ref. [3]. Multinucleon 
mechanisms (not properly QE) in the neutrino nucleus scattering were studied in Ref. |18j] . The approach is based on 
a systematic many body expansion of the W boson absorption modes that includes one, two and even three nucleon 
mechanisms, as well as the excitation of A isobars followed by their absorption and the production of pions. The 
parameters of the calculation had been previously fixed in the study of photon, electron and pion interactions with 
nuclei. The model successfully describes jl^ the published MiniBooNE neutrino CCQE-like data, finding a large 
contribution of two nucleon excitation mechanisms or, using the usual many body naming, 2 particle - 2 hole (2p2h) 
channels. In a latter work 20], we investigated the energy reconstruction procedure with the aim of determining 
the influence of the multinucleon excitations on the extraction of neutrino energy unfolded cross sections from the 
measured flux-average data. We found that the 2p2h contribution strongly distorts the shape of the total CCQE-like 
flux unfolded cross section for neutrino beams. 

In this work, we apply the same model (loj for QE plus multinucleon mechanisms to antineutrino nucleus scattering 
and compare with the recent MiniBooNE data [2^ . The differences with the neutrino calculation are straightforward. 
The sign of the antisymmetric piece of the leptonic tensor is changed and that implies a change on the sign of the parity 
violating pieces. In terms of the nuclear structure functions, this means a change of sign in the pieces proportional to 
W3 in the differential cross section, following the notation of Ref. [5|] . The boson selfenergy in the nuclear medium 
is replaced by that of the . This is achieved by exchanging the role of protons and neutrons in all formulae. Finally, 
proper care has to be taken with the energy balance because of the different Q values. We will always use the same 
axial mass of our previous papers {Ma — 1.049 GeV). 

The results presented in this paper correspond to QE and 2p2h contributions to the f-^ C — > /i+ + X process. 
Scattering on other nuclei like ^^O or ^^Fe produce similar cross sections per nucleon. Apart from the D, v energy 
dependence of the cross sections and their ratio and the muon angular distribution, we will also show results of the 
double differential cross section (Pa /dcosO^dTf^ folded with the MiniBooNE flux. 




E (GeV) 

FIG. 1: Cross section contributions divided by the number of protons for v -|-^^ C scattering as a function of the antineutrino 
energy. Curves labeled appx correspond to the approximated unfolded cross section for the MiniBooNE flux. Experimental 
points from Ref. [2^ . 
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The results are shown in Figs. 1-4. In Fig. [1] we show the energy dependence of the QE and the 2p2h contributions 
to the Dfj^ +^'^ C cross section. The cross sections are smaller than for neutrinos. In this case, 2p2h channels also 
provide an important contribution which cannot be neglected. One of the relevant questions is to assess the quality of 
the standard procedures of energy reconstruction. Our results are similar to those obtained for the neutrino case. In 
the figure, we have labeled appx the approximated unfolded cross sections^. We observe that for genuine QE events 
the curve is hardly modified. However, for 2p2h events there is a substantial displacement towards low energies that 
is also reflected in the summed cross section. In other words, some "high" energy antineutrinos would be identified as 
"low" energy ones. Although this statement is of a general validity, the detailed distortion of the cross section curve 
depends on the antineutrino flux. The full model for the unfolded aappx agrees well with the experimental data ^, 
whereas the pure QE contribution predicts lower values for the cross section. 
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FIG. 2: Ratios of neutrino versus antineutrino cross sections in ^^C as a function of the energy. Solid line: QE + 2p2h. 
Dash-dotted line: QE. Triangles: Ratio obtained with the MiniBooNE flux using the unfolding procedure of Ref. [201 ]. 

Fig. [21 shows the ratio between the neutrino and the antineutrino induced reactions. The inclusion of the 2p2h 
mechanisms slightly reduces the ratio (by less than a ten percent at the peak) without much affecting the shape. Thus, 
in our model the relative importance of the 2p2h channel is somehow larger for antineutrinos. A similar trend, although 
with a stronger reduction, has been found by Amaro et al. [isj in the superscaling approximation and by Meucci et 
al. [m in the EDAI option of the relativistic Green's function model, which agrees well with the i^fj, MiniBooNE data. 
Other works like Ref. Tg^, which reaches agreement with MiniBooNE QE neutrino data by modifying the magnetic 
form factors of the bound nucleons, and Ref. lead to larger values for this ratio. This latter reference is of particular 



^ See Ref. [20| | for details on the unfolding procedure 

^ For this observable, there is an additional 17.2% normalization uncertainty not shown in the figure. 
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interest as it is the closest to our approach and contains a microscopical evaluation of 2p2h mechanisms similar to 
ours. However, there are some important differences which amount to a more comprehensive inclusion of mechanisms 
in our scheme and some approximations used in the calculations. A detailed discussion on these differences can be 
found in Ref. Finally, we find that the unfolding procedure respects the Oyjoi, ratio for energies above 0.6 GeV 
but produces a strong distortion at lower energies. 

In Fig.[3l we present the muon angular distribution for 1 GcV neutrinos (antineutrinos). The v scattering is more 
forward peaked than the v one for both the QE and the 2p2h contributions. This is due to the different sign of the 
piece proportional to the W3 structure function (cf. Eq. (10) of Ref. Q) which leads to a negative interference for 
the V case. 
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FIG. 3: Muon angular distributions for 1 GeV neutrinos (antineutrinos) scattering on ^^C. 



As discussed in Ref. [19|, excluded the observation of final hadrons, a good observable to compare with theoretical 
models is the double differential cross section corresponding to the energy and angle of the final muon, because it 
avoids the problems of energy reconstruction and depends just on directly measured quantities. In Fig. 21 we show 
our results for that cross section for the muonic antineutrino flux of MiniBooNE. In general terms, we successfully 
describe these data. For all angular windows 2p2h mechanisms play an important role and lead to a much better 
agreement with data. Their contribution tends to concentrate at low muon energies so that more energy is transferred 
to the nucleus than in pure QE events. Comparing this plot with its neutrino counterpart [sl. [l9l|. we observe a quite 
different pattern because of the more forward peaked character of the antineutrino induced reaction. There are also 
some changes concerning the antimuon energy distribution that are mostly due to the lower average energies of the 
MiniBooNE v flux. 

In summary, we have studied non mesonic antineutrino nucleus cross sections in a model that starts from a relativistic 
local Fermi gas description of the nucleus, includes RPA correlations and multinucleon effects finding a good agreement 
with experimental data. The same model had been quite successful for the neutrino cross section and contains no 
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FIG. 4: Muon angle and energy distribution a / dcosd ^dT^ per proton on a ^'^C target folded with the MiniBooNE muon 
antineutrino flux. Different panels correspond to the various angular bins labeled by their cosinus central value. Black-solid 
line is the full model including QE (relativistic and with RPA) and 2p2h mechanisms. Red-dash-dotted curve corresponds to 
QE and blue-dashed curve to 2p2h events. Data from Ref. [231, errors only account for the shape uncertainty. Besides, there 
is an additional normalization uncertainty of 17.2%. 



free parameters. The effects of multinucleon events are seen to be larger, relative to the QE contribution, than for 
the neutrino case and considerably improve the agreement with experiment. We have also found that the size of the 
2p2h cross section limits the reliability of the unfolding procedures based on the assumption of QE kinematics for the 
non mesonic events. 
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